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Abstract 

Vector-like quark is a common feature in many new physics models. We study the vector- like 
quark production as an s-channel resonance at the Large Hadron Collider for the vector-like quarks 
mainly mixing with the top- and bottom-quark in the Standard Model. We emphasize that the 
leptonic angular distribution can be used to discriminate various vector- like quark models. 
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I. INTRODUCTION 



Vector-like quark (vector-quark) exists in many models of new physics (NP) beyond 
the Standard Model (SM), e.g. extra dimensional models, Little Higgs models, dynamical 
models, etc. The vector-quark, denoted by Q, could mix with the SM quarks through 
Yukawa interaction. In particular, the mixing of vector-quark with the third generation 
quarks in the SM is quite common in the Little Higgs models [l| and composite Higgs 
boson models {2]. Several possibilities for their electroweak quantum numbers are listed 
in Ref. jsj] (see Table HJ. The phenomenology of the vector-quark production in hadron 
collisions has been studied either in the pair production 141 or in the association production 

n n 

with a SM quark |5|, [6[. In the work we consider the mono vector-quark production as 
an s-channel resonance via strong magnetic q-g-Q couplings at the Large Hadron Collider 
(LHC). Depending on whether the vector-quark would mix with SM chiral quarks, we can 
further focus on the decay channel that the vector-quark decays into a top-quark plus a 
jet, i.e. Q — > tg and other electroweak decay channels. The top-quark polarization can 
be measured through the charged lepton angular distribution from top-quark decay inside 
the top-quark's rest frame and the polarization is sensitive to the chirality of the g-t-Q 
coupling. The SU(2) x U(l)y charged and neutral currents will be modified according to the 
quantum numbers of vector-quarks after adding Yukawa mixings, therefore the branch ratios 
for vector-quark decay to electroweak gauge bosons vary in several vector-quark models. 
Discovering the vector-quark signal and further measuring the top-quark polarization would 
help us to distinguish several new physics models. 

The q-g-Q coupling is forbidden by the Ward Identity at the tree level. The interaction 
could be induced by the new physics (NP) at high energy scale (A). Rather than working 
in a ultra-violet completion theory, we adapt the approach of effective field theory (EFT) in 
this study. The vector-quark might carry different quantum number under the electroweak 
symmetry of the SM (57/(2) l x U (l)y). It could be a weak isospin singlet, doublet or triplet. 
Below we list all the possible gauge invariant dimension-6 effective operators which describe 
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TABLE I: Electroweak quantum numbers for vector-quark multiplets, which could mix with the 
SM quarks through the Yukawa interaction. 



SU(2)t 



Y 



Vn 



V x V Y T x T Y 



U D 




X 

u 

D 







u 

D 
Y 



2/3 -1/3 1/6 7/6 -5/6 2/3 -1/3 



the strong magnetic g-q-Q interaction: 

K u,R 



C 



gqQ 



— q L a^\ A S u ^G A v + -jf q L a^\ A S d ^G A v 



+ 
+ 



K>X,L 



V 2L a^\ A u R 4>G 



A 

/IV 



Kd,L 
A 2 



V 2L a^\ A d R cj ) G 



A 

fW 



— V XL ^X A u R <pG A u -r 

KX,R 



V YL a^\ A d R (j)G 

K Y,R 



/us 



+ ^jf Txro^\ a (0rV) G% + ^ T Y RO^\\^q L )G 



A 2 



(1) 



where (j) is the SM Higgs boson doublet, G A is the field tensor of gluon and A is the cut 
off scale of our effective theory. Note that the r 1 denotes the weak isospin in the basis 
of (+,0,—). Table |T] shows the detailed gauge quantum number for all the vector-quarks 
scenarios. 

The paper is organized as follows. In Sec. II we consider the simplest scenario where the 
vector-quark interacts with light quark only through the dimension-6 operators in the sense 
that it won't mix with the standard model chiral quarks. The top-quark is predominately 
left-handed polarized when it originates from a gauge singlet or triplet vector-quark. On the 
other hand, it is mainly right-handed polarized when it comes from a gauge doublet vector- 
quark decay. Since top-quark is the only "bare" quark in the SM, we consider the decay 
mode of Q — > tg to utilize the top-quark polarization to distinguish the gauge singlet /triplet 
and gauge doublet models. In the section III , we add the Yukawa interaction to couple the 
vector-quark to the SM chiral quarks such that fermions with the same quantum number 
can mix. Electroweak precision measurement i.e. the T parameter and Z — > bb, are adopted 
to constrain the parameter space. After that, we proceed to analyze the channel of pp — > 
Q — > bl + u. We derive the leptonic angular distribution with respect to the gluon moving 



FIG. 1: Feynman diagram of the process of ug — > Q — > gt(—> be + v e ). 

direction in the center of mass frame of initial partons, which can be used to discriminate the 
singlet/triplet model from doublet model under certain assumptions. We are going to show 
that the mixing pattern for each vector quark scenario determines their discovery potentials 
in the LHC. 

II. EXCITED QUARKS AND COLLIDER SIMULATION 

We are interested in the following process that the mono produced vector-quark decays 
into a top-quark plus a jet ( see Fig. CD for Feynman diagram ): 

ug -> Q->tg, t-> b£ + u, (2) 

with Q carrying an electromagnetic charge 2/3. In this section we assume that Q does not 
mix with the chiral top quark in the SM such that they can only decay through dimension-6 
operators. For clarity, we only include the positively charged state in the following discussion. 
Furthermore, the leptonic decay mode of the top-quark is considered in the analysis because 
the top-quark spin is maximally correlated with the charged lepton £ + . The drawback is the 
neutrino is invisible. It escapes the detection and yields a signature of large missing energy 
($t) & t the collider. One has to determine the neutrino momentum to fully reconstruct the 
top-quark kinematics, which is the key for the top-quark polarization measurement. 

After the spontaneously breaking of the electroweak symmetry, the effective operators 
described by Eq. [T] induces an effective chromo-magnetic-dipole coupling of g-u-Q&s follows: 

L D UlPl + f R P R ) QG, u {k) + h.c. (3) 

where v = 246 GeV denotes the vacuum expectation value for the Higgs field and Pl, Pr 
are the chiral projectors. We further assume that the heavy quark Q interacts with all three 
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up-type quark in the SM (u, c, t) universally. As we can see from its original dimension-6 
form, for each vector-quark scenario the interaction should be either left handed or right 
handed depending on their representations in the SU(2) gauge group, i.e. fi and can 
not be present simultaneously. This property makes the dimension-6 operators distinct 
from the dimension-5 operators written down in the excited quark model |7| • Similar FCNC 
anomalous magnetic dimension-6 operators with the heavy quark Q replaced by the top 
quark are adopted to study the single top production js-10], whose coupling bounds are 



recently explored in the ATLAS experiments 11 ). 



Since the excited heavy quark Q has an electric charge 2/3, we denote it as one heavy 
vector like T-quark in the following. The single T-quark production cross section at the 
LHC is : 

°{qg ->■ T) = - / dxidx 2 {f g /p(xi)f g /p(x 2 ) + (xi <-> x 2 )} 5(s - m 2 T )-^r(fl + fn), (4) 

where fi/p is the so-called parton distribution function (PDF) describing the possibility of 
finding a parton % inside a proton with a momentum fraction Xj. In the simplified model, 
the decay width is : 

r(T "> qg) = 6^A%4 > (5) 

which clearly displays that the branching ratio of T — > tg increases with while the 
branch ratios of T — >■ ug , eg decrease with rrir, and they all approach to 1/3 in the large 
limit. Since our effective operators contribute to the quark-gluon final states in the hadron 
collider, the upper limit of the dijet cross section put constraints on the couplings J'l and 
Jr 12, uM- The cross section a times the CMS acceptance factor A for pp — > T — > ug, eg 
at a 7 TeV LHC is compared with the corresponding CMS upper limit for the quark-gluon 
final state as shown in Fig. [2) When we choose the cut-off A to be 1.2 TeV, the magnetic 
coupling strength need to satisfy the conditions of fi < 3.0, fn = (or f L = 0, /r < 3.0). 
In this paper, we consider the relatively small magnetic coupling scenario. The inclusive 
cross section for the process of qg — > T — > bC + vg with £ + = e + , /i + when the center of mass 
energy is 7 TeV, 8 TeV and 14 TeV is also plotted in Fig. (2) We choose the couplings to be 
f L = 0.5, f R = (or f L = 0, f R = 0.5) with the cut-off A fixed to be 1.2 TeV. 

The top-quark polarization can be best measured in its leptonic decay mode t — > M + u. 
The collider signature of interest to us is one charged lepton, two jets plus large missing 
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FIG. 2: The upper panel plot is the inclusive cross section (in the unit of picobarn) times CMS 
acceptance A for the process pp — > T — > qg at a 7 TeV LHC. The red dotted line is for ug and eg 
final states via one excited T-quark, the blue solid line is for tg final state via one excited T-quark 
and the dark red line is the CMS upper limit for quark-gluon final states. The couplings are chosen 
as fi = 3.0 and Jr = while the cut-off scale A is fixed at 1.2 TeV. The lower panel plot is the 
inclusive cross section for the process pp — > T — > qg at 7 TeV, 8 TeV and 14 TeV LHC in the green 
solid line, red dotted line and blue solid line respectively. The couplings are chosen as Jl = 0.5 
and Jr = with the same cut-off scale A = 1.2 TeV. 

energy, where the large missing energy originates from the invisible neutrino and one of 
the two jets is the gluon in association with top-quark production and the other one is the 
b-quark from top-quark decay. The signature suffers from a few SM backgrounds as follows: 
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Single-t production: single top-quark can be produced via the electroweak interaction 
in the Standard Model. It proceeds through the s-channel decay of a virtual W 
(qq 1 — >■ W* — > tb), the t-channel exchange of a virtual W boson (bq — > tq' , bq' — > tq), 
and the associated production of a top-quark with two jets (qq ; , gq —> tjj). Note that 
tjj includes the process of gb — > tW~ and with W~ subsequently decay into two jets. 
The single top plus one jet process is the intrinsic backgrounds as it yields exactly 
the same signature as the signal event. On the the hand, the top quark associated 
production with two jets process is the non-intrinsic background, since this process 
can only mimic the signal in the case that one of the two jets produced in association 
with the top-quark escapes the calorimeter detection. 

W+ two jets production processes: among them are the Wbj production and the Wbb 
production which have the same signature as the signal events; the other one is Wjj 
production process but it requires one of the two jets to fake the fe-jet. As to be shown 
later, it still contributes as the largest background even after including the small faking 
efficiency Additional small background is from WZ production with Z gauge bosons 
decay into two jets. 

ti production: it is a non-intrinsic background because, in order to mimic the signal, 
either the two jets from the i hadronic decay or the charged lepton from the i leptonic 
decay need to get lost in the detector. 



Both the signals and backgrounds are generated by MadGraph/MadEvent |14J. The 



Q. 



CTEQ6L1 parton distribution functions [15| are used in this study. When generating the 
Wjj backgrounds, we impose soft cuts on the transverse momentum of the both jets, i.e. 
PtU) > 5 GeV, to avoid the collinear singularity. Note that such soft cuts do not affect our 
analysis as we impose much harder cuts on both jets in the following analysis. The inclusive 
cross section for the signal event i.e. pp — >• T — >• tg with the following decays t — > bt + v 
(£ + = e + , is shown in the second column in Table [TXJ For illustration we choose six 

benchmark masses for the T quark and choose fi = 0.5 and = 0. The cut off scale is set 
to be A = 1.2 TeV throughout this work. Other results of different couplings can be easily 
obtained by rescaling: 

<r,(h, fn) = <Ts(h = 0.5, f R = 0) 4 (f L + f 2 R ). (6) 
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TABLE II: Cross section (in the unit of fb) of the signal events for various tut at the 14 TeV 
LHC (upper panel) and of the SM backgrounds (lower panel). The model parameters are chosen 
as f L = 0.5 , f R = and A = 1.2 TeV. 



m T (GeV) 


no cuts 
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AM W < 10 GeV 
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500 GeV 


1893.51 


512.192 


256.287 


240.759 


233.659 


600 GeV 


1453.86 


417.329 


322.756 


304.365 


295.569 


700 GeV 


1120. 


335.439 


292.542 


277.59 


268.687 


800 GeV 


869.751 


276.362 


255.054 


242.617 


235.703 


900 GeV 


680.2 


223.344 


211.95 


200.387 


194.299 


1000 GeV 


540.21 


181.403 


174.866 


166.736 


161.658 




Backgrounds 


no cuts 


AM W < 10 GeV 

basic cuts Pf > 200 GeV H T > 400 GeV 












AM t < 10 GeV 


Wjj 


1.8532 * 10 7 13476.5 


637.501 


633.794 


163.082 


Wbb 


48054. 


1191.74 


13.215 


12.014 


3.604 


Wbj 


2774.2 


36.758 


0.277 


0.277 


0.139 


WZ 


544.13 


99.494 


0.449 


0.422 


0.027 


tj 


29619. 


4214.78 


311.0 


294.71 


281.38 


tjj 


28234. 


5522.57 


364.219 


290.81 


237.166 


tb 


1025.1 


247.649 


12.945 


12.004 


9.718 


tbj 


27745 


975.237 


40.230 


31.907 


26.358 


ti(b£u) 


25431. 


118.254 


3.178 


1.272 


0.0 


total 


1.86954 * 10 7 25883 


1383.01 


1277.21 


721.474 



To simulate a realistic detection, we impose the event-selection cuts as follows: 

p b T > 20 GeV, | jfc | < 2.5, 
Vt > 20 GeV, 1 77, | < 2.5, 

ARjj > 0.4, ARji > 0.4 , (7) 
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where AR(= a/ (Ar/) 2 + (A</>) 2 ) is the separation between any two observable final state 
particles (not including neutrinos) , and A0 and At] are the separation in azimuthal angle 
and rapidity, respectively. In order to let those non-intrinsic background to mimic the signals, 
additional veto cuts are demanded for jets or leptons which are not detected (either falling 
into a large rapidity region or carrying a too small transverse momentum to be detected), 

veto cuts: PtU,^) < 10 GeV or \q{j,^)\ > 3.5 . (8) 

with t] denoting the rapidity of the jets or leptons in the final state. The vetoing cut is 
imposed at the same time with the basic cut when we are performing the events selection 
and the corresponding cross section after the event selection is presented in the third column 
of Table [Til We model the detector resolution effects by smearing the final state energy 
according to 

-FT = 7 , © B, 9 
E y/E/GeV 

where we take A = 10(50)% and B = 0.7(3)% for leptons (jets). In addition, we require 
that one jet is tagged as 6-jet with a tagging efficiency is 50%. We also apply a mistagging 
rate for charm-quark e c ^ = 10% for Pt(c) > 50 GeV. The mistagging for a light jet is 
e u ,d,c,s,g^b = 0.67% for p T (j) < 100 GeV and 2% for p T (j) > 250 GeV. For 100 GeV < 
Pt{j) < 250 GeV, we linearly interpolate the fake rate given above. 

At this stage of the analysis, the background rate is one or two order of magnitude larger 
than the signal rate. Moreover, the dominate background comes from the Wjj process, 
followed by the single-top plus jets production and other small background processes. In 
order to study the efficient cuts that can significantly suppress the background rates while 
keeping most of the signal rates, we examine the px distributions of the leading-p T jet as 
shown in Fig. EJ The leading jet in the signal originates from the light jet produced in 
association with the top-quark. Owing to the large mx, the leading jet peaks in the large pt 
region. On the contrary, the leading jets in the backgrounds, either from the QCD radiation 
or from the VF-boson decay, tend to peak in the small pt region. Such a distinct difference 
enables us to impose a hard cut on the first leading jet, 

PrUist) > 200 GeV, (10) 

to suppress the huge backgrounds. In the fourth column of Table |TTJ we show the cross 
sections after the above cuts. This cut increases the signal-to-background ratio by a factor 
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FIG. 3: The left panel is the px distribution of the leading jet in the signal and background events 
after the basic and veto selection cuts. The right panel is the Ht distribution of all the final objects 
in signal and background events after the basic and veto selection cuts. Each distribution plot is 
normalized. 

of 17.2 while keeping about 92% of the signals for my = 800GeV. The biggest reduction 
in the background rate comes from the Wjj production, but all the other backgrounds are 
reduced sizably as well. 

One can further impose a hard cut on the Ht variable, the scalar sum of the transverse 
momentum of all the measurable objects in the final states i.e. Ht = Y2Pt+ with % 
sum over the visible objects. The Ht distribution of the signal events peaks in the large 
Ht > 400 GeV region, while the one of the SM backgrounds is mostly located in the small 
H T < 400 GeV region (see Fig [3]). But this cut does not influence both the signal and 
background too much after the hard Pt(Jut) > 200 GeV cut . 

We are interested in measuring the top-quark polarization, which requires a full recon- 
struction of the top-quark kinematics. One confronts the invisible neutrino in the final state. 
Assuming the missing transverse momentum {$t) comes entirely from the neutrino, i.e. 

p v (x) = - $ T (x), Pu(y) = - #t(j/), (11) 

the longitudinal momentum of the neutrino can be reconstructed by the on-shell condition 
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FIG. 4: Contour of the discovery potential of the signal at the 14 TeV LHC with an integrated 
luminosity of 1 fb _1 . Both the 3 a exclusion limit S/y/S + B = 3, (red-dashed line) and the 5 a 
discovery limit S/\/B = 5 (blue-solid line) are shown in the figure. 

of VF-boson: 

m 2 w = {p e +p,) 2 . (12) 

The quadratic equation yields a two-fold solution. We only pick the real solutions and 
abandon the complex ones due to the W gauge boson's width effect. Since the VF-boson 
comes from a top quark, both solutions are used to reconstruct the top quark mass, so 
that we can pick the one which gives a mass closer to 173 GeV. The neutrino momentum 
reconstruction is conducted after the large px cut. After neutrino reconstruction we are able 
to determine the kinematics of both the PF-boson and top quark. Two mass-window cuts 
around mw and m t are imposed to optimize the signal events, 

AM W = \m lv - m w \ < 10 GeV, AM t = \m M - m t \ < 10 GeV. (13) 

Fig. |4] plots discovery potential of the signal event in the plane of tut and fi {fn = 0) 
after the mass-window cuts. We display both the 5cr discovery limit and the 3a exclusion 
limit in the significance plot, which shows that it is promising to observe the single T-quark 
production at the LHC with an integrated luminosity of 1 fb _1 . 

The behavior of the angle between the lepton in the rest frame of the top quark and the 
top quark moving direction in the center of mass frame can be factorized out and described 
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FIG. 5: The cos 9 distribution of the signal with no cut (dashed line) and after the mass-window 
cuts (solid line): (a) the singlet/triplet T-quark, (b) the doublet T-quark . 



by the following equation [lg, ll 

1 da 1 ,„ _ . ,„ 

- = - l + acos0 £ , 14 

a a cos Bg 2 

where a = 1 for the pure left-handed heavy top quark and a = — 1 for the pure right-handed 
heavy top quark. Figure |5] displays the cos 6g distribution with no cut and after imposing 
the mass- window cut. 



III. VECTOR-QUARK AND SM QUARK MIXING 

In the previous discussion we assume the vector-quarks do not mix with the SM quark at 
the tree level. However one could write down a gauge-invariant Yukawa interaction between 
the vector-quark and the SM quarks as follows: 

£y = -y u qLH c u R - y d q L Hd R - X u q L H c U R - \ u V 2 lH c u r - \ d V 2 LHd R 
-X u V XL Hu R - \ u q L r a H c T£ R - \ u q L T a HTJ R - MU L U R 
-MV 2L V 2R - MV XL V XR - MTxlTxr - MT Y lT yr , (15) 

The Yukawa interaction generates a mixing between the SM quarks and the vector-quarks 
at the tree level after the spontaneous symmetry-breaking. The singlet vector-quark and the 
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triplet vector-quark exhibit a similar mixing pattern, while the doublet vector-quark have 
a different mixing pattern. For simplicity, we consider the scenario that the vector-quarks 
mix only with the third-generation quark in the SM . Consider one pair of vector-quark T' L R 
or one pair of vector-quark B' L R which will mix with the chiral top-quark t' L R or the chiral 
bottom quark b' L R in the following way : 

t' L , R = cosei R t LjR + smei R T L}R} (16) 
T' LjR = - sin 6l R t L , R + cos 6l R T LiR , (17) 

b' LR = cos6 d L>R b L , R + sm6 d LtR B L!R , (18) 
B' L)R = - sin 9 d LtR b LiR + cos 6l R B LtR , (19) 

where t^R and Tl ,r labels the physical top-quark and heavy vector-quark, respectively. 
Defining x t = \ u v/ \/2 and Xf, = Xdv/^/2, the mixing angle can be calculated by diagonalizing 
the mass matrix. For the singlet /triplet vector-quark we obtain 

3in^ (d) = — , (20) 

Bin^ (d) = — ■ (21) 
while for the doublet vector-quark 

sm^) = — ( 22 ) 

sm^ = "7^ ==■ (23) 

Since the mixing would inevitably modify the W-t-b and Z-b-b couplings in the SM, the 
electroweak precision measurements at the low energy would severely constrain the mixing 



Mx t{b) 


\A M2 


- m\ b) f + M^x\ h) 




m t (b)X t (b) 


aA m2 


-m 2 m f + MW m 




m t (b)Xt(b) 


\A M2 


-ml {b) f + M*x\ b) 




MXtQb) 


\A M2 


-m 2 m ) 2 + M^ t{b) 



parameters [181 ] - Before proceeding with the further analysis, we are going to consider the 
low energy precision test in several vector-quark models and put bounds on the parameter 
space using current experimental constraints. When fixing X},, the parameter describing the 
mixing of the bottom-quark and the vector-quark, the stringent constraint in the (M, x t ) 
parameter space of the singlet and doublet vector-quark models is from the T-parameter, 
while it is from the Z-b-b coupling for the triplet vector-quark model. The mixing of vector- 
quark with standard model chiral quark modifies the charged current of W± gauge bosons 
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and the neutral current of W% gauge bosons simultaneously and the mass splitting among 
the fermions will contribute to T parameter through the vacuum polarization of the gauge 
bosons 
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20] . For the singlet vector-quark U\, two types of doublet vector-quark T> 2 



and T>x and the triplet vector-quark 7x, their specific contributions to T parameter are 
formulated in Eq. 

3 



AT; 



IQits 2 w c 2 w 



[sin 2 fie+iyr, y b ) - sin 2 9^9 + (y t , y b ) - cos 2 9 L U sin 2 %9 + (y T , y t )} (24) 



[sin 2 eZe+iyr, y b ) - sin 2 9 L J+{y u y b ) 



+ (sin 2 9 L U + sin 2 9*)9 + (y u y x ) + (cos 2 6 L U + cos 2 9*)9 + (y T , y x ) 

+ 2 sin 9 L U sin 6*6- (y t , yx) + 2 cos 9 L U cos (y T , y x ) 

- (4 cos 2 ^ sin 2 9 L U + cos 2 ^ sin 2 6*)6 + {y u y T ) 

- 4 cos ^ sin ^ cos ^ sin ^(y*,^) 



(25) 



r> 2 



1671-%^ 



cos 



sin 



sin 



2 „2 



[(cos 2 (^ - ^) + sin 2 9* sin 2 0* - l)9 + (y t , y b ) 



9 L d ) + cos 2 0* cos 2 ef)e + (yT,y B ) 
9 L d )+ cos 2 0* sin 2 9«)9 + (y t ,y B ) 
9 L d )+ sin 2 0* cos 2 9«)9 + (y T ,y b ) 
+ 2 cos(0£ - 0^) (sin 0* sin 9*9- (y t , y b ) + cos 9* cos 0^0_ (y r , y fl )) 
+ 2sin(^ - d L )(cos0«sin0*0^, y B ) - sin 9* cos 9*9 _(y T ,y b )) 
■ cos 2 9* sin 2 0*0+(yt, 2/t) - cos 2 9* sin 2 0^ + (y 6 , y B )] 



(26) 
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AT Tx 



((cos 9 L d cos % + V2 sin 9 L d sin + 2 sin 2 Q* sin 2 0* - l)0 + (y t , y b ) 

+ 2^(005 ^ cos 9% + v 7 ^ sin 0^ sin 0£) sin 0^ sin 0f 0_ (^, y 6 ) 

+ ((cos 0^ sin 0^ - v 7 ^ cos 9 L d sin 0^) 2 + 2 cos 2 0^ sin 2 9*)9 + (y t ,y B ) 

- 2^(005 0^ sin 9 L d - V2 cos 0^ sin 6%) cos 0^ sin 0f 0_ (y t , y B ) 

+ ((- V2 cos 0^ sin 0^ + cos 9 L d sin 0^) 2 + 2 cos 2 0* sin 2 9*)9 + (y T , y b ) 

- 2y/2(- \Pl cos 0£ sin 0^ + cos 9 L d sin 0^) cos 9* sin 9%9_{y T , y b ) 

+ ((sin ^ sin ^ + y/2 cos 0, L cos 0^) 2 + 2 cos 2 0* cos 2 9*)9 + (y T , y B ) 

+ 2v / 2(sin 0^ sin 0^ + v 7 ^ cos 9 L d cos 0£) cos 0f cos 9*9- (y T , y B ) 

+ (2 sin 2 ^ + 2 sin 2 9*)9 + (y t , y x ) + 4 sin 9 L U sin 0^0_ (y t , y x ) 

+ (2 cos 2 ^ + 2 cos 2 0^)0+ (y r , y x ) + 4 cos ^ cos 9*9 L (y T , y x ) 

- (cos 2 9 L d sin 2 d L + 4 cos 2 0* sin 2 9 R )9 + (y b , y B ) 

- 4cos0£sin0£cos0>n0*0_(2/&,2/ B ) -cos^ 2 sin^ 2 + (2/r, 2/t)] (27) 

where 2/; = mf/m 2 z is the dimensionless fermion mass squared rescaled by the inverse mass 
squared of the Z gauge bosons. I have already checked that all the divergence are can- 
celled due to the unitary property of mixing matrix and the relation between those mass 
eigenstates. Since the T parameter measures the effect of custodial symmetry breaking, 
unctions 0+(2/i, y 2 ) and y 2 ) are zeroes when yi = y 2 and they are defined by the 

20| : 
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9 + {y 1 ,y 2 ) = ?/i + 2/2 log — 

Vi ~ 2/2 yi 

at \ o / — fyi + yz, Vi 

M2/i,2/2) = 2 V2/i2/2 log— - 

\V\-V2 2/2 



(28) 
(29) 



Another source of T parameter contribution comes from the Higgs mass deviation from the 
its reference value mh re f = 120 GeV: 



AT, 



167TC 2 



log 



in 7 



m 



(30) 
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The electroweak precision measurement gives the fit for the T parameter is T = 0.05 ± 
0.11 23[. Assuming the Higgs mass is 125 GeV 2lJ, |22[ , and adding two sources of 
T parameter contribution, the contour of T-parameter in (M, x t ) parameter space for the 
singlet vector-quark and doublet vector-quarks are plotted in Figure Eta-d). For the singlet 
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FIG. 6: (a) T parameter constraint for M and xt in the singlet IA\ model ; (b) T parameter 
constraint for M and x± in the doublet T>x model; (c) T parameter and i?6 constraints for M and 
X[, in the doublet T>2 model with xt fixed to be 80, 90, 100, 120 GeV ; (d) T parameter constraint 
for M and x t in the doublet T>2 model with x\, fixed to be 20 GeV. 

scenario as depicted by Fig. Efa), only the Higgs gives a very small negative contribution and 
the singlet top quark can achieve larger positive contributions, therefore we get a positive 
T-parameter deviation in the interested vector-quark mass region. Requiring T < 0.17 and 
varying the gauge invariant mass M in the range of (500 — 1000) GeV, the upper limit for 
the top quark mixing parameter x t should be in the range of (97.4 — 153.4)GeV. Since the 
mixing pattern and hyper charge assignment determine the SU (2) gauge bosons currents, we 
see another situation in the T>x doublet vector quark scenario. The Fig. E(b) illustrates the 
T parameter constraint for the T>x doublet model with the same mixing parameter x t as the 
singlet model. The extra 5/3 charged heavy quark X in the doublet model does not mix with 
any standard model quarks, but they will contribute to the AT© X . The heavy quark loop 
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FIG. 7: T parameter and Rb constraints for M and xt in the triplet 73c vector-quark model . 

effects can achieve notable positive T parameter deviation when x t is large while they give 
notable negative T parameter deviation when the x t is small. Requiring —0.07 < AT < 0.17, 
two stringent green bands in the (M, Xt) plane are the permitted parameter space. We are 
more interested in the lower green band located in the small Xt region. As the gauge invariant 
mass varies in the range of (500 — 1000) GeV, the negative T parameter bound constrains 
the Xt to be in the range of (54 — 78) GeV. For the doublet model with one heavy top 
quark and one heavy bottom quark, there will be two independent mixing parameters Xt 
and Xb for the up type quark and down type quark respectively. The T parameter mainly 
constrains the x t and the other mixing parameter Xb is constrained by the z — > bb which will 
be considered following. In the Fig. [6](c), the cyan curve, green curve, red curve and the 
blue curve depict the upper bound T = 0.17 in the (M, Xb) parameter space for four sample 
values Xt = 80, 90, 100, 120 GeV . As we can see, in the low Xb region which is allowed by 
the Rb constraints, the T parameter is not sensitive to the Xb but sensitive to the other 
mixing parameter Xt and the upper bound for gauge invariant mass M will increase as we 
enhance the value of x t . Fig. [6]^d) plots the ATx> 2 dependence on M and x t with the fixed 
value Xb = 20 GeV. Similar to the singlet model the T parameter mainly receives positive 
contribution from the heavy fermions in the T> 2 doublet model. After imposing the upper 
bound of T < 0.17, we can find that Xt is constrained to be in the range of (81 — 123) GeV 
in correspondence with the M in the range of (500 — 1000) GeV. Finally, we discuss the 
situation in the 73c triplet model. Only the heavy top quark and the heavy bottom quark 
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in the triplet model will mix with the SM chiral quarks and the two mixing parameters are 
related to each other by virtue of Xb = v2a; t . The T parameter dependence for M and Xt is 
plotted in Fig. which shows that it almost does not put any constraint in the parameter 
space compared with the tight Rb requirement. 

In the situation with one pair of heavy bottom quark present, we need to consider the 
constraints from the Z — > bb. Both the corrections to Z-b^-bi couplings and the corrections 
to Z-bR-bR couplings in the doublet model as well as in the triplet model are induced through 
the Yukawa mixing between the vector bottom quarks and the chiral bottom quarks. The 
Sg^ p and 5g^ p can be translated to be the deviation 5Rb by the following formula Eq. [3T] : 

5R b = 2R b (l - R b ) (9*^ + _9R 5g sp\ (31) 

9l = — (-5 + ^s w ) g R = ^—s w , (32) 

where and gn are Z gauge bosons couplings to the left handed bottom quark and the right 
handed bottom quark in the Standard Model. Rb is defined as Y{Z — > bb)/T(Z — > hadrons) 
with its SM value Rb = 0.21578^q;qqq8. The deviation SRb due to the new physics effects is 
constrained by electroweak measurements to be {23 1 : 

5R b = 0.00051 ±0.00066 (33) 

We can write down the Z-b^-b^ couplings and the Z-b^-bR couplings in the T> 2 doublet model 
and in the l~x triplet model. For the doublet model we have: 

9zbL = —{-\ + \s 2 w)=9T b L (34) 

9zm = -\~ sin 2 0* + = sin 2 # + g ThR (35) 

while for triplet model we have : 

9zbL = -\^^e L d +M-\ + \sl / ) = - l -^^e L d+ gr bL (36) 

— 9 „: T .2/ifl 1 1 9 2 _ 9 , 2 nR , n sm /o 7 \ 

9zbR sm V d + - — s w sm V d + g ZbR (37) 

Substituting Eq.[34l- 137] into Eq. [31], we get the Rb deviation due to the vector-quarks 
in the doublet model and in the triplet model respectively, which in turn gives stringent 
constraints for Xb and x t in each scenario. As we can see from Fig. [6](c) which depicts the 
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FIG. 8: Feynman diagram for single-T production and decay (ug — > T — > bW + — > b£ + V£). 

situation for the T>2 doublet model, the negative Rb bound limits the x& parameter to be in 
the range of (23 - 46) GeV as M varies from 500 GeV to 1000 GeV. Fig. [7] shows that it 
is the upper limit of Rb < 0.00117 which put constraints on M and x t for the l~x triplet 
model. The green band is the allowed parameter space, i.e. for 500 GeV < M < 1000 GeV, 
we need 20 GeV < x t < 39 GeV. 

In this section we still consider the T-quark is produced via anomalous q-T-g couplings, 
but it could decay through the Yukawa interaction shown above. Note that after the mixing 
the T-quark as an s-channel resonance is no longer purely polarized. After adding the 
Yukawa mixing, the T-quark gets three additional decay modes opened: th, tZ and bW + . 
Here we proceed to analyze the process of gq — > T — > bW + with a subsequent PF-boson 
leptonic decay (see Fig. [8] for Feynman diagrams), yielding a collider signature of b£ + $t- 
The main SM backgrounds are (1) W + j where the light jet can mimic bottom quark; (2) 
W + bj as well as t-channel single-top production qb — > q't — > W + bj, in both cases we require 
that one of non-b tagged jets escapes detecting; (3) s-channel single-top production channel, 
qq' — > tb — >• W + bb, where one of the bottom quarks escapes the detection. The basic cuts 
to trigger the events are: 

p\ > 50 GeV, \r] b \ < 2.0, 

p\ > 20 GeV, 1 7ft | < 2.4, (38) 
AR bl > 0.7. 

The cut table for the signal event (singlet case and doublet case) at the bench mark point 
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TABLE III: Cross section (fb) of the signal process pp — > T — > bW + — > bl + v and SM backgrounds 
for various kinematics cuts. The top table is for the singlet vector-quark with the choice of /x, = 0, 
ffi = 0.2 and tut = 800 GeV. The middle table is for the doublet vector-quark with the choice of 
/l = 0.3, /r = and tut = 800 GeV and with an addition input Xb = 20 GeV. The bottom table 
is for the SM backgrounds. 



xt (GeV) no cuts basic cuts P\ 


> 200 GeV ARui < 1.5 A 


m w < lOGeV 


40 GeV 211.8 


73.643 


69.841 


66.283 


66.283 


60 GeV 223.1 


77.059 


73.121 


69.596 


69.596 


80 GeV 228.02 


80.046 


76.056 


72.237 


72.237 


100 GeV 230.44 


80.216 


76.518 


73.199 


73.199 


120 GeV 232.18 


81.890 


77.595 


74.019 


74.019 




xt (GeV) no cuts b 


iasic cuts Pt 


> 200 GeV AR ul < 1.5 A 


,mw < lOGeV 


40 GeV 166.25 


57.847 


54.921 


52.336 


52.336 


60 GeV 114.36 


39.694 


38.088 


36.115 


36.115 


80 GeV 88.792 


30.314 


28.968 


27.552 


27.552 


100 GeV 75.728 


25.369 


24.426 


23.211 


23.211 


120 GeV 68.41 


22.599 


21.70 


20.732 


20.732 




SM bg no cuts 


basic cuts P t 6 > 200 GeV AR ul < 1.5 . 


Am w < lOGeV 


w + j 1.95 * 10 7 


5889.06 


175.97 


175.97 


175.97 


w+bj 2774.2 


80.32 


1.942 


1.248 


1.248 


tj 29619 


5386.22 


42.947 


7.405 


7.405 


tb 1025.1 


22.14 


0.28 


0.22 


0.22 


total 1.95334 * 10 7 


11377.7 


221.139 


184.843 


184.843 



in? = 800 GeV as well as for the background is shown in Table IIHI The third column shows 
the numbers of event after the basic cut. We demand 6-tagged when imposing the basic cut 
where the veto cut for the non-intrinsic background events are also included. The 6-tagging 
requirement effectively reduces the Wj background by a factor of 0.5 x 10~ 3 , while it still 
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keeps about 1/3 signal events. In the signal event the 6-jet comes from the heavy vector- 




FIG. 9: (a) The pt distribution of the leading jet after the basic and veto selection cut ; (b) the 
distribution of the A.R(l + v) between the lepton and the reconstructed neutrino after the basic and 
veto selection cut. Each distribution is normalized. 

quark decay. It carries a large px- On the contrary, the fe-jets in the background events, 
either from the gluon splitting or from the top-quark decay, exhibit a relatively soft px- 
The difference can be seen from Fig. |9]^a) where the pt distribution of the 6-jet is plotted. 
Following the basic cut and 6-tagging, we impose a hard px cut on the 6-jet, 

Pr(b) > 200 GeV. (39) 

As shown in Table II III the cut reduces the background more than one order of magnitude, 
but it leaves the signal events almost untouched. The reason that we do not consider the 
non-intrinsic background events from W + bb and W + Z is that bb from either gluon splitting 
or Z bosons decay can not simultaneously pass the veto cut as well as the hard px{b) cut . 

To fully reconstruct the signal event, the kinematics of the invisible neutrino is to be 
determined from the VF-boson on-shell condition. We pick up the solution with a smaller 
magnitude in the two-fold solutions. The H^-bosons from the heavy vector-quark T decay 
is highly boosted such that its decay products, the lepton and neutrino, tend to collimate 
and yield a small AR(£ + , v) separation as we can see from the blue-solid curve in Fig[9](b). 
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On the other hand, the background events are more evenly distributed in the relative large 
AR(£ + , v) region. 



Significance : S/V^ Significance: S/Vb" 




500 600 700 800 900 1000 500 600 700 800 900 1000 



M (GeV) M (GeV) 

FIG. 10: Significance contour in the (x,M) parameter space: (a) the singlet IA\ model (fn = 0.2); 
(b) the doublet V 2 model (f L = 0.3 and x b = 20 GeV); (c) the doublet V x model (f L = 0.6) ; (d) 
the triplet Tx model (/# = 0.3). 

Fig. [10] displays the significance contours at the 14 TeV LHC for several vector-quark 
models. The doublet T>y is ignored since there is no heavy T-quark in that model, and 
the triplet Ty models is not considered either because it has a very small branching ratio 
for the decay of T — > W + b . We will emphasize some characteristic in those contours and 
illustrate the reasons as follows. The Ui, T>x and Tx models exhibit a similar pattern of 
the discovery potential, i.e. xt increases with M for a fixed significance. While the doublet 
T>2 model shows a different pattern as Xt decreases with M. The difference is caused by 
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the decay branching ratio of the vector-quark in various models [18j. The branching ratio 
of T — > W + b decay always increases as the value of Xt is enhanced in the U\, T>x and Tx 
models. However in the T> 2 model the corresponding branching ratio can be approximated 
by x\J{x\ + 1x\) in the large M limit. There is a tension exists between Xj, and x t and 
the branching ratio of T — > W + b decreases with x t for a fixed Xf, value, which leads to the 
slopping pattern displayed in Fig. [TQT b). As we can see from the significance plot for the 
U\ singlet model with ff> = 0.2, that the significance contour is not sensitive to the mixing 
parameter in the large Xt region, such that the a = 5 discovery limit curve clearly cuts the 
M to be less than 800 GeV . Note that the T>x doublet model has a smaller discovery 
potential as indicated by its magnetic coupling. The reason is that the branching ratio of 
T — > W + b is suppressed by a factor of m 2 /M 2 in that model as compared with the other 
three ones. 




FIG. 11: The cos0£ distribution between the lepton and gluon moving directions with no cut 
(dashed line) and after the mass-window cut (solid line): (a) the singlet T-quark, (b) the doublet 
T-quark (when there is no heavy bottom quark or its effect can be ignored). 

It is convenient to analyze the charged-lepton angular distribution between the lepton and 
the gluon moving directions in the center of mass frame, which can be used to distinguish 
the chiral property of T-b-W couplings. The differential cross sections with respect to cos Qg 
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for the ug — > b£ + u in vector-quark models in the limit of mw <C ^ are found to be: 
1 da(u R g -j> bt+v) g^ 2 g* 2 

-z ~ n = -ry- ^ (1 - COS 6» f ) + — rr - i jj-j (1 + COS0*) 

a dcos0 e 9^ +9w 9w +9w 

R 2 

0(m 2 w /s)- cos 9 e (40) 



9y/ 2 9w 2 



1 da(u L g b£+u) g^ 2 g% 2 

- = -tt~, — Wli 1 + cos °e) + -ft; — r^(I-cos^) 

a dcosOt 9w +9w 9w + 9w 

R 2 

+ lIZ. R 2 -0(m 2 w /s).cos9 e (41) 
9w i" S'vk 

Note that the Oim^j s) correction is only for the term proportional to the right-handed 
T-b-W coupling. The detail derivation of exact results for the leptonic angular distribution 
in the rest frame of the heavy top quark is put in the Appendix. Fig. [TT] displays the cos 0$ 
distribution with no cut and after the mass-window cut. For the singlet U\ model, since 
the anomalous u-g-T interaction is right handed and the coupling of T-b-W is purely left 
handed, the angular distribution of cos6^ should be proportional to (1 — cos6')/2 as indicated 
by Eq. (40] ; On the other hand, the anomalous u-g-T interaction is left handed and the 
coupling of T-b-W is also purely left handed in the doublet T>x model, therefore the lepton 
angular distribution should be proportional to (1 + cos6 l )/2, which is in coincidence with 
the analytic result (see Eq. [41]). The situation is more complicated in the doublet V 2 
model. In that model we have g^/9w ~ i x t m t) I i^bM) such that the angular distribution 
shows the superposition of (1 =F cos6 l )/2. When Xb ^> x t , the coupling of T-b-W is mainly 
right-handed and the angular distribution is similar to the singlet U\. In the other limit of 
x t ^> Xb, the decay of T — > bW + is possible to be dominated by the left-handed coupling as 
long as M is not too large, so that its angular distribution is similar to the doublet T>x- The 
leptonic angular distribution should serve as an effective analyzing power when the V — A 
chiral structure of the T-b-W vertex is not too much modified. 



IV. CONCLUSION 



In this paper, the single heavy top quark production via strong magnetic interaction 
is studied and we explore the possibility for vector-quarks to be discovered in both the 
anomalous decay and electroweak decay channel. The leptonic angular distribution is a 
favored analyzing power for identifying the top quark polarization and distinguish varieties 
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of vector-quark models. We use the collider simulation to explain that the top polarization 
in the channel of T — > tg is determined by the chirality property of the excited quarks. The 
Yukawa mixing does not change the situation since the mixing only dilutes the corresponding 
branch ratio in each specific model. However the V — A chiral structure of T-b-W is possible 
to be modified by those Yukawa mixing . We conduct a detail analysis for the parameter 
space constrained by the electroweak precision test in various vector-quark models. The 
channel of T — > bW + has very less standard model backgrounds, and the leptonic angular 
distribution in that channel can at least be utilized to distinguish the U\ singlet model from 
the T>x doublet model. Once we found the signals of heavy T vector-quark and measured 
its chirality, we are capable to reconstruct its mass as a resonance in a single production 
process. Although we have shown that the LHC is sensitive to the signals of vector like 
quarks in the single production channel, the cross section of such process depends on the 
coupling strength of strong magnetic interactions, otherwise the pair production of heavy 
top quarks should be a better discovery channel . 
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Appendix A: Leptonic angular distribution in cm. frame 



We derive the moving direction of charged lepton relative to the gluon moving direction 
in the center of mass rest frame for the process ug — > T — > b£ + u . The z axis is defined by the 
moving direction of the charged lepton £ + in the cm. frame of the initial partons . Since the 
production and the decaying can be factorized using the narrow width approximation, we 
further have the mass on shell relation i.e. k^y = M^, such that the energy and momentum of 
the bottom quark is determined by energy-momentum conservation to be: Eyy = + 
and p = ^(1 — ^f L ). The four momentum for the initial partons and the bottom quark, 
the lepton and the intermediate W + gauge bosons can be expressed as: 

P 9 = (^,-^sin^,0,^cos^) (Al) 

Pu = ^ sin 0,, 0,-^ cos 0,) (A2) 

k b = (j},psm9 bi cos4> M ,psmO M sm(f) bi ,pcos9 b i) (A3) 

k w = (E w , -p sin 9 be cos (p b£ , -p sin 9 M sm(p M ,-p cos 9 b£ ) (A4) 

k e = (E e ,0,0,E e ) (A5) 

where 9 £ is polar angle for the incoming gluon makes with the lepton moving direction (z- 
axis), 9 b i is angle between the bottom quark and the lepton and the 4>m is the corresponding 
azimuthal angle. In the left-handed and right-handed strong magnetic interaction cases, the 
amplitude squared is calculated to be: 

64 

\M(u L g be + v)\ 2 = -j^glv 2 f 2 L p g ■ p u (^ 2 s k e ■ p u k ve ■ k b + gy\ 2 Ml k ve ■ p g k e ■ k b ) 

nS(k 2 w - m 2 w ) C F ■ N c 

T w m w ((s-M 2 ) 2 + T 2 T M 2 ) 1 ' 

64 

\M(u R g -> be + v)\ 2 = -^gjv 2 f^p g ■ p u (g^M^ K ■ p g k ve ■ k b + g% 2 s k ve ■ p u k e ■ k b ) 

ixb{k 2 w - m 2 w ) C F ■ N c 

T w m w {{s-M 2 Y + T 2 T M 2 ) 1 ) 

where Cf = 4/3 and N c = 3 are the color factors. A symmetry exists between the fi and 

fn scenarios that k ve is exchanged with k e and the (V-A) coupling g^ is exchanged with the 

(V+A) coupling gy\r. Phrasing the final states phase space in the rest frame of the heavy 

top quark, the differential cross section is expressed as: 

1 IMP 4 

da = — jTgTg j^XE dE b dE e d cos M<M06i (A8) 
26 



with the following integration limits: 



2 2 



< E b < 
-1 < cos8 f < 1, < (i 



< 2tt 



(A9) 



We have used the narrow width approximation for the W + gauge bosons when we are 
calculating the amplitude squared and the delta function can be written as : 



5{k 



2 

W 



m 



w) 



2^ \ 2 



(A10) 



The cos &u is not independent which can be written in terms of E b and E £ by the following 
expression: 

s-2yfi(E b + E £ ) : 1 



cos 9, 



(All) 



2E b Ei 

Now we put all the elements into the amplitude squared and after integrating out E b , E# 
and two azimuthal angles of <pe and <p b £ we get the following differential cross sections with 
respect to cos 0£ which depict the angular distribution of the charged lepton in the center of 
mass rest frame: 



da{u L g — > bt+v) 
d cos Qe 



1 1 



3 2 2 n 7r 2 A 4 Y w m w s 



- g^ 2 M 2 ( (m 2 w - s) 2 {2m 2 w + s) + Y2m% / s - 12m 2 w s 2 + 12m^slog 



w 



+ g& 2 M 2 (m 2 w - s) 2 (2m 2 w + s] 



((s - M 2 ) 2 + YD 



cos9e 
(A12) 



da(u R g -> M+v) g 2 2 v 2 f 2 R 



dcos9 e 3 2 2 n 7r 2 A 4 T w m w s 

, 2 



i (g^ 2 ■ M 2 (m 2 w - s) 2 (2m 2 w + s) 



+ gw 2 s[ {m 2 w - s) {2m 2 w + s) + 12m^s - \2m 2 w s 2 + 12m^s log 

1 



in 



9w 2s ( m w 



,)'(2 



m 



ir 



1) 



1 — cos 6e) 
cos ^ 
(A13) 



(( s -M 2 ) 2 + r 

As we can see, only the term proportional to the (V-A) coupling has exactly the 1 ± cos 9g 
distribution and the term proportional to the (V+A) coupling is slightly corrected to be 
(1 =F cosQe) ± Oijnlyl ' s) cos^ , the deviation can be ignored in the limit of m w <C y/s . 
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